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Oh, We show that enhanced extra mixing in low-mass red giants can result in a 

fluorine abundance that is correlated with abundance variations of other elements 
"^ [ participating in H burning, such as C, N, O and Na. This finding is used to 

explain the fluorine abundance variations recently found in bright red giants of 

^ , the globular cluster M4. 
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Introduction 



o 

H i In many galactic globular clusters, there are star-to-star abundance variations of C, 



N, O, Na, Mg and Al, while the heavier elements show almost no dispersion (see reviews 
by Denissenkov 2004b; Gratton et al. 2004, and references therein). The signs of these 
variations, their pairwise correlations, and the constancy of the sum of C+N+0 in some 
^ ' clusters (Pilachowski 1988; Smith et al. 1996; Ivans et al. 1999) indicate that we are most 

likely seeing the by-products of H burning in the CNO-, NeNa-, and MgAl-cycles. Re- 
garding their production place and process, the most plausible alternatives proposed so far 
include the so-called hot-bottom burning of hydrogen at the bottom of convective envelope 
in intermediate-mass (4 < M/Mq < 6) asymptotic giant branch (IM-AGB) stars (D'Antona 
et al. 1983; Ventura et al. 2001), and the hydrogen shell burning in low-mass (M < 2 Mq) 
red giant branch (RGB) stars (Sweigart & Mengel 1979; Denisenkov & Denisenkova 1990; 
Langer & Hoffman 1995). The former case would imply primordial abundance variations 
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in the RGB stars where the abundance anomahes are observed, while the latter explana- 
tion would require in situ mixing. The physical interpretation of these results is therefore 
significant for understanding the origin of these anomalies. 

Recently, Smith et al. (2005a) have added ^^F to the list of nuclides whose abundances 
vary from star to star in globular clusters. They have found a deficit of the fluorine abundance 
([F/Fe]^ < 0) anti-correlated with the abundance of Na in 7 bright red giants of the globular 
cluster M4. The range of the fluorine abundance variation is ~ 0.8 dex. Since ^^F is effectively 
destroyed in the hot-bottom H burning in the reaction ^^F(p,Q;)^^0, while low-mass stars have 
so far been considered as its net producers (Forestini et al. 1992; Jorissen et al. 1992), Smith 
et al. (2005a) have interpreted their measurements as an indication that the star-to-star 
abundance variations in globular clusters originated in the IM-AGB stars. 

In this paper, we give an alternative interpretation of the new observational data by 
demonstrating that low-mass RGB stars could produce the F-Na anti-correlation in exactly 
the same way as they contribute to the global 0-Na anti-correlation (Denissenkov & Van- 
denBerg 2003; Denissenkov et al. 2006). 



2. IM-AGB and RGB Pollution Scenarios 

In some globular clusters, the 0-Na and Mg-Al anti-correlation and the N-Na correla- 
tion have been found not only in evolved RGB stars but also in main-sequence (MS) turn-off 
and early subgiant stars. Because these stars have interior temperatures that are too low for 
the NeNa- and MgAl-cycles to work, they must have abundance variations that predate the 
RGB phase; these are either present ab initio or reflect accretion of material from IM-AGB 
stars (Briley et al. 1996; Gratton et al. 2001; Ramirez & Cohen 2002; Cohen et al. 2002; 
Briley et al. 2002; Grundahl et al. 2002; Harbeck et al. 2003). 

In IM-AGB stars, the envelope material processed in hot-bottom H burning is immedi- 
ately delivered to the interstellar medium via a strong stellar wind. On the other hand, in 
RGB stars the H-burning shell is separated from the bottom of the convective envelope by a 
radiative zone. However, extensive spectroscopic data on the surface chemical composition 
of low-mass RGB stars as a function of luminosity, supported by appropriate stellar models, 
show that most of these stars experience some extra mixing that connects the H burning 
shell with the convective envelope (Sweigart & Mengel 1979; Charbonnel 1994, 1995; Gratton 



^We use the standard spectroscopic notation: [A/B] = log (iV(A)/iV(B))-log (A^(A)/iV(B))Q, where N{A) 
and A^(B) are number densities of the nuchdcs A and B. 
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et al. 2000; Denissenkov & VandenBerg 2003). Near the RGB tip, before undergoing the 
core He-flash and becoming a horizontal branch star, a low-mass red giant loses a consider- 
able amount of its envelope mass with a stellar wind (D'Cruz et al. 1996), thus depositing 
products of the H-shell burning to the interstellar medium. In globular clusters, both the 
IM-AGB stars and the low-mass stars more massive than the present-day MS turn-off stars 
(i.e. 0.9 M0 ^ M < 2 Mq) have already completed their lives. So, either of them could 
pollute the globular-cluster interstellar medium with the ashes of H-burning. 

Unfortunately, neither the AGB nor the RGB pollution scenario can quantitatively 
reproduce all of the observed star-to-star abundance variations in globular clusters simulta- 
neously. For example, when destroying ^^O in the IM-AGB stars the hot-bottom burning 
depletes ^^Mg even more (Denissenkov & Herwig 2003). It also keeps [C/Fe] > —0.5 (Denis- 
senkov & Weiss 2004). Besides, the third dredge-up in AGB stars should increase the sum 
of C+N+0 (Fenner et al. 2004). None of these theoretical predictions is supported by ob- 
servations. At present, IM-AGB star models are being modifled towards complying with 
constraints imposed by the element abundance anomalies in globular clusters (e.g., Ventura 
& D'Antona 2005). 

In low-mass RGB stars, extra mixing starts to manifest itself when the H-burning shell, 
advancing in mass, erases the chemical composition discontinuity left behind by the bottom of 
convective envelope at the end of the first dredge-up (Gratton et al. 2000; Shetrone 2003). At 
this moment, the evolution of red giants slows down for a while, which produces a prominent 
feature (the RGB bump) in the differential luminosity functions of globular clusters (Zoccali 
et al. 1999; Riello et al. 2003). Therefore, extra mixing in low-mass stars is said to work on 
the upper RGB, above this bump luminosity. Below the bump luminosity, on the lower RGB, 
extra mixing is thought to be shielded from the H-burning shell by a gradient in the mean 
molecular weight associated with the composition discontinuity (Sweigart & Mengel 1979; 
Charbonnel, Brown & Wallerstein 1998; Denissenkov & VandenBerg 2003) or to operate very 
slowly (Chaname et al. 2005; Palacios et al. 2006). 

The main problem of the RGB pollution scenario is that for the majority of upper 
RGB stars, the observed pattern of surface abundance anomalies can only be produced by 
in situ mixing that only penetrates to the outer part of the H-shell, where the CN branch 
of the CNO-cycle is operating. This is what Denissenkov & VandenBerg (2003) have called 
"canonical extra mixing". They demonstrated that its depth and rate (diffusion coefficient) 
can be parameterized by any pair of correlated values within the close limits specified by 
AhgT ^ 0.19 and D^ix ~ 4 x 10^ cm^ s-\ to AhgT ^ 0.22 and D^^^ f^ 8 x 10^ cm^ s^\ 
Here, A log T is the difference between the logarithms of temperature at the base of the 
H-burning shell and at the maximum depth of extra mixing. Thus, canonical extra mixing 
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can be responsible for the evolutionary decline of [C/Fe] in upper RGB stars, both in the 
field and in star clusters. It can also explain the decrease in surface Li and ^He abundances 
with increased L, the strong reduction of the ^^C/^^C isotopic ratio and the slight increase 
in the N abundance, but it does not affect O, Na, Mg and Al. 

However, if extra mixing in upper RGB stars penetrates the H-burning shell deeper 
than in the canonical case, it could dredge up material deficient in O and enriched in Na 
(Denisenkov & Denisenkova 1990) and even in Al, under certain assumptions (Langer & 
Hoffman 1995; Cavallo et al. 1996; Denissenkov & Tout 2000; Denissenkov & Weiss 2000). 
Denissenkov & VandenBerg (2003) have proposed that in some upper RGB stars canonical 
extra mixing may be switched to its enhanced mode with much faster and somewhat deeper 
mixing. If the extra mixing in these stars is driven by differential rotation of their radiative 
zones then such enhanced extra mixing could be caused by their spinning up as a result of 
tidal synchronization in close binaries (Denissenkov et al. 2006). 

There is also an unexplored possibility that the depth and rate of canonical extra mixing 
do not remain constant along the whole upper RGB but increase toward its tip. This 
hypothesis is supported by the following arguments. Firstly, observations show that in some 
globular clusters the anti-correlated abundance variations of C and N in red giants become 
larger when stars approach the RGB tip. Moreover, the values of [N/Fe] > 1 in some of 
these stars indicate the dredge-up of material in which not only C but also a fraction of O 
has been converted into N (Smith et al. 2005b). Secondly, at least in the globular cluster 
M13, the relative number of upper RGB stars with the 0-Na anti-correlation increases 
with luminosity (Johnson et al. 2005). Note that all of the metal-poor field stars used by 
Denissenkov & VandenBerg (2003) to calibrate the parameters of canonical extra mixing have 
log L/ Lq < 2.8, i.e. they are located a magnitude below the RGB tip. Hence, if canonical 
extra mixing got enhanced within the last magnitude of the upper RGB then Denissenkov 
& VandenBerg (2003) would not have seen it. 



3. ^^F and Extra Mixing in Upper RGB Stars 

3.1. Model Parameters for the M4 Red Giants 

The stellar evolution code, input physics and simple diffusion model of extra mixing in 
upper RGB stars used here have been described by Denissenkov et al. (2006). The code has 
since been developed to allow studies of the pre-MS and horizontal branch (HB) evolution. 
As in many other similar codes, our zero-age HB models are constructed using information 
about the internal structure of the RGB tip models in which the core He-flash just sets in. 
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We have found that for a model star with the initial mass M = 0.83 M©, helium and 
heavy-element mass fractions Y = 0.24 and Z = 0.002 our RGB evolutionary track, zero-age 
HB and HB evolutionary tracks fit the color-magnitude diagram (CMD) of evolved stars in 
M4 reasonably well (Fig. 1). Our adjusted stellar parameters give a theoretical metallicity 
[Fe/H] Ri — 1.0 and an age of evolved stars ~ 14 Gyr. In Fig. 1, we have applied the cluster 
reddening E{B — V) = 0.40 and the distance modulus (m — M)v,o = 12.83. Our adopted 
parameters are in good agreement with those used by Alcaino et al. (1997) who fitted the 
M4 main-sequence CMD with a 14 Gyr isochrone of Bergbusch & VandenBerg (1992) for 
Y = 0.2388, [Fe/H] = -1.03, E{B -V) = 0.42 and (m - M)y,o = 12.80. The spectroscopic 
metalhcity of M4 is [Fe/H] = -1.18 (Ivans et al. 1999). 



3.2. ^F and Canonical Extra Mixing 

Abundance profiles of some nuclides participating in the H-shell burning, including 
^^F, in a bump luminosity model star are plotted in Fig. 2. The vertical line marked 
A log T = 0.19 is shown at a depth characteristic of canonical extra mixing. The nearly 
vertical line at the left of the plot shows the profile for H; this coincides with the lower edge 
of the H-burning shell. Canonical extra mixing occurs when the H-burning shell erases the 
composition discontinuity at the mass coordinate M,. ^ 0.266 M©. Most or all stars more 
luminous than this would show a reduced abundance of ^^C (Denissenkov & VandenBerg 
2003). However, we see that canonical extra mixing will not change the surface abundance 
of i^F. 



3.3. ^^F and Enhanced Extra Mixing 

Note that all of the 7 red giants studied by Smith et al. (2005a) are located within 
one magnitude of the RGB tip {asterisks in Fig. 1). As we have already mentioned in §2, 
there are some arguments supporting the idea that at least the depth of canonical extra 
mixing increases in the vicinity of the RGB tip. In particular, the decline of [C/Fe] becomes 
stronger within the last magnitude of the RGB, which is accompanied by extremely high N 
abundances ([N/Fe] > 1) supposedly signifying the dredge- up of material with O partially 
processed into N (Smith et al. 2005b). Smith et al. (2005a) claim that they have found 
a similar decline of the C abundance A(^^C) = log A^(C)/A^(H) + 12 with the absolute 
bolometric magnitude in their M4 stars {circles with errorbars in the left panel in Fig. 3). 

In order to test the hypothesis that canonical extra mixing gets enhanced near the RGB 
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tip, we have proceeded as follows. First, starting with a model at \ogL/LQ = 1.75 (this 
value is slightly less than the bump luminosity; it corresponds to Mboi = 0.37, My = 0.72, 
and V = 13.55), we have computed its evolution up to logL/L© = 2.74 (Mboi = —2.10, 
My = —1.37, and V = 11.46). In this computation, the depth and rate of extra mixing had 
their canonical values A log T = 0.19 (Fig. 2) and -Dmix = 4 x 10^ cm^ s~^, or logDmix = 8.6. 
Abundance profiles in the radiative zone in the final model are shown in Fig. 4. After that, 
we have continued our stellar evolution computations toward the RGB tip with increased 
depth and rate of extra mixing. The parameter A log T has been reduced to 0.05, while for 
logDmix we have considered the following higher values: 9.2, 9.3, 9.4, 9.5, and 9.6. Already 
looking at Fig. 4, one can conclude that such enhanced extra mixing will dredge up material 
in which deficits of ^^O and ^^F should correlate with overabundances of N and Na. Curves 
in left panel in Fig. 3 demonstrate how the surface C abundance declines with Mboi in our 
models with enhanced extra mixing. In Fig. 5, a combined evolution of [C/Fe] is shown as 
a function of My for the post-first-dredge- up {dotted line), canonical extra mixing {dashed 
curve), and enhanced extra mixing {solid curve) phases. Here, points with errorbars are 
data for bright red giants of the globular cluster NGC 7006 from the work of Smith et al. 
(2005a), the 4 brightest of them having [N/Fe] > 1.3 (given the large uncertainties of the 
[C/Fe] data in NGC 7006, the difference in metallicities [Fe/H]]\^QQ7QQg— [Fe/H]]y[4 ^ —0.5 
is unimportant for the comparison). 

Interestingly, if we take the fluorine abundances measured by Smith et al. (2005a) in the 
M4 red giants at their face values and plot them as a function of Mboi, we find a correlation 
of A{^^F) with Mboi {circles with errorbars in right panel in Fig. 3) that looks none the 
worse than the correlation between A{^'^C) and Mboi (left panel). Even more interestingly is 
that the curves in the right panel representing results of our evolutionary computations with 
enhanced extra mixing could nicely reproduce that correlation if it were real. We also find 
a theoretical anti-correlation between [Na/Fe] and [F/Fe] that is qualitatively similar to the 
observational one revealed by Smith et al. (2005a) (Fig. 6). 

Therefore, alternatively to the IM-AGB pollution scenario, we propose that the fluorine 
abundance variations in the red giants of the globular cluster M4 may actually be a signa- 
ture of enhanced extra mixing in low-mass RGB stars. Moreover, we strongly encourage 
spectroscopists to increase the size of the sample of bright red giants with known fluorine 
abundances in globular clusters because if A{^^F) indeed correlates with Mboi this will be a 
direct evidence of enhanced extra mixing in these stars. On the other hand, if the fluorine 
abundance is not found to decline with increasing luminosity near the RGB tip this will not 
necessarily reject our hypothesis that its deficit is due to enhanced extra mixing in upper 
RGB stars. There will still remain a possibility that the star-to-star abundance variations 
of ^^F, like those making up the 0-Na anti-correlation, were produced by enhanced extra 
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mixing in the low-mass red giants that had been shghtly more massive than the present-day 
MS turn-off stars in globular clusters and that had completed their lives in the past (the 
RGB pollution scenario). 

We have also examined the possibility that the F/Na anticorrelation could have been 
produced by larger-than-expected temperature errors. The F abundances were computed 
using the temperatures determined by Ivans et al. (1999), who obtained them by measuring 
the line depth ratios of temperature-sensitive species following Gray & Johanson (1991). 
The stated random errors in the temperature were about 50 K. Using the sensitivity of the 
abundances to temperature, we find that the dispersion in the F or Na abundances could 
have been produced if the random errors in temperature were as large as 125 K. Since both 
lines have the same sense of temperature dependence, however, such large random errors 
would produce an erroneous correlation between the F and Na abundances, instead of the 
observed anticorrelation. 

Note that for a Salpeter initial mass function the total mass lost by those red giants 
before they arrived at the zero-age HB is comparable to the mass delivered to the interstellar 
medium by the IM-AGB stars (Denissenkov 2004a). This estimate assumes that every upper 
RGB star loses 0.2 M© before it undergoes the core He-flash. For comparison, our HB 
model stars with {B — V) < 0.7 on the CMD of the globular cluster M4 in Fig. 1 all have 
M < O.6IM0, i.e. each of them has lost more than 0.22 Mq. It is far from clear, however, that 
the problem of depositing the lost mass from either AGB or RGB stars onto the other stars 
in the system after they reached the MS has been satisfactorily addressed. For this reason, 
determining the contribution of direct mixing on the observed abundances is important. 



4. Conclusion 

In this work, we have shown that the reduced abundances of ^^F (i.e. [F/Fe] < 0) anti- 
correlated with [Na/Fe] found by Smith et al. (2005a) in the red giants of M4 do not rule 
out the RGB pollution scenario in which canonical extra mixing switches to its enhanced 
mode. The enhanced extra mixing could result from tidal spin-up of upper RGB stars in 
close binaries (Denissenkov et al. 2006). Another possibility is that canonical extra mixing 
gets enhanced toward the RGB tip due to some internal physical processes in single stars. 
We emphasize that there is observational evidence of this, one of which could be a correlation 
of A(^^F) with Mboi like that plotted in the right panel in Fig. 3, provided it is confirmed 
in future spectroscopic observations. At present the reduced fluorine abundance cannot be 
considered as a strong argument in favor of the IM-AGB pollution scenario as the only 
plausible interpretation of star-to-star abundance variations in globular clusters. On the 



contrary, if A{^^F) is correlated with Mboi in low-mass stars near the RGB tip, this will 
support the original hypothesis of Denisenkov & Denisenkova (1990) that extra mixing in 
upper RGB stars can penetrate the H burning shell deep enough to dredge-up material 
enriched not only in N but also in Na, and deficient not only in C but also in O and F. 
A small number of stars in the fast evolutionary phase, at the RGB tip, could explain the 
absence of field upper RGB stars in which the anti-correlated O and Na surface abundances 
are produced and dredged up in situ (Gratton et al. 2000). 

We acknowledge support from NASA grant #NNG05 GG20G. 
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Fig. 1. — CMD of evolved stars in M4 from the data of Cudworth & Rees (1990) [open 
circles), Smith & Briley (2005) {filled circles) and Smith et ah (2005a) (asterisks). Curves 
represent our evolutionary tracks for a model star with M = 0.83 M©, Z = 0.002 and 
Y = 0.24 computed from the pre-MS to the RGB tip and, after that, from the zero-age HE 
(starting with M = 0.58 Mq and M = 0.62 Mq) to the beginning of the early AGB phase. 
The cluster reddening E{B — V) = 0.40 and the distance modulus [m — M)v,o = 12.83 have 
been applied. These values are close to those used by Alcaino et al. (1997) and others. 
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Fig. 3. — The correlation between the surface C abundance and Mboi claimed by Smith et al. 
(2005a) to exist in their observed bright red giants of the globular cluster M4 {circles with 
errorbars in the left panel) and a possible correlation between the fluorine abundance and 
Mboi revealed by us in their published data (the right panel). Curves represent results of 
our evolutionary computations with enhanced extra mixing near the RGB tip for the mixing 
depth Alog T = 0.05 (Fig. 4) and rates log Dmb. = 9.6, 9.5, 9.4, 9.3, and 9.2 (from the 
bottom to the top curve). 
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Fig. 4. — The same as in Fig. 2 but for a model near the RGB tip. For enhanced extra 
mixing with the depth shown with the sohd hne, the surface fluorine abundance undergoes 
considerable changes (see the right panel in Fig. 3). 
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Fig. 5. — Evolutionary decline of [C/Fe] with My caused by canonical extra mixing with 
A log T = 0.19 and log -Dmix = 8.6 {dashed curve) followed by enhanced extra mixing with 
A log T = 0.05 and log -Dmix = 9.6 {solid curve). Dotted line gives the post-first-dredge-up 
value of [C/Fe]. Circles with errorbars are observational data for bright red giants of the 
globular cluster NGC7006 (Smith et al. 2005a). 
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Fig. 6. — The Na-F anti-correlation revealed by Smith et al. (2005a) in 7 bright red giants 
of the globular cluster M4 {circles with errorbars). Curve shows our theoretical dependence 
of [F/Fe] on [Na/Fe] obtained as a result of enhanced extra mixing with A log T = 0.05 and 
log -Dmix = 9.6 that is assumed to replace canonical exra mixing when stars approach the 
RGB tip (for details, see text). 



